Prictional coupling between sliding and spinning motion 
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We show that the friction force and torque, acting at a dry contact of two objects moving and 
rotating relative to each other, are inherently coupled. As a simple test system, a sliding and 
spinning disk on a horizontal flat surface is considered. We calculate, and also measure, how the 
disk is slowing down, and find that it always stops its sliding and spinning motion at the same 
moment. We discuss the impact of this coupling between friction force and torque on the physics 
of granular materials. 



Sliding friction and incomplete normal restitution are 
normally the main dissipation mechanisms at the contact 
between two solid grains. They are largely responsible 
for the fact that the flow properties of granular media 
differ from those for liquids and solids. Their microscopic 
origins are currently under intense investigation (see e.g. 
jjj ) . On large scales compared to the grain diameter they 
sometimes transform into unexpected phenomenological 
friction laws, which have recently been discussed ^ [j| |J, 
H ||. As the influence of rolling and torsion friction is 
commonly regarded as negligible, they have been much 
less investigated so far. However, it turns out that in 
certain situations they may become crucial, for instance 
for the stabilization of pores in cohesive powders j?], [|. 
Another striking phenomenon will be discussed below. 

In fact, rolling and torsion friction are indispensable 
for a unified view of the dissipation mechanisms at the 
contact of two viscoelastic spheres, because on the one 
hand incomplete restitution and rolling friction, and on 
the other hand sliding and torsion friction are coupled, 
which will be the main point of this letter. These four 
dissipation mechanisms correspond to the six degrees of 
freedom of the relative motion at the contact between two 
solid spheres. The relative motion of two solid spheres 
has three translational degrees of freedom, characterized 
by a velocity vector with one normal component v n (de- 
formation mode) and two tangential components v t (slid- 
ing mode), and three rotational ones, characterized by an 
angular velocity vector, again with two tangential com- 
ponents u>t (rolling mode) and one normal component 
u n (spinning mode). While the viscoelastic dissipation 
mechanism of normal restitution and rolling friction cou- 
ples v n and uj t || |l^| , the dissipation due to sliding and 
torsion friction couples vt and to n . 

In this Letter we focus on the coupling between v t 
and uj n . For viscoelastic spheres these are the domi- 
nant dissipation channels in the quasistatic limit, where 
(1 — restitution coefficient) [jll], p"2|, as well as the co- 
efficient of rolling friction vanish 19|. The reason, why 
torsion friction, i.e., the torque leading to a decrease of 
u! n , is often neglected, is that it involves the radius of the 
contact area between the two spheres and hence is small. 



FIG. f : (a) A sliding and spinning disk on a flat horizontal 
surface, (b)-(d) The relative velocity field on the surface of 
the disk at e = 0.2, e = 1, and e = 5, respectively (e = v/Rui). 



Therefore, in order to make our point more clear, instead 
of the contact between two spheres we consider a flat disk 
on a horizontal flat surface with nonzero initial transla- 
tional and angular velocity. The disk is lying on one of 
its sides, and we assume that this side is in full contact 
with the table during the motion [see Fig. ga)] . The fric- 
tion force and torque acting on the disk will slow down 
the sliding and spinning motion until the disk stops mov- 
ing. We address two questions: (i) how are the friction 
force and torque related to each other, (ii) what does this 
imply for the coupling of sliding and spinning motion? 

First we calculate the friction force and torque acting 
on the disk as a function of its instantaneous velocity 
and angular velocity. We apply the Coulomb friction 
law, which says that the magnitude of the friction force 
is proportional to the normal force, while its direction 
is the opposite of the direction of the surfaces' relative 
velocity. Assuming that the mass distribution of the disk 
is homogeneous, the friction force is 
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where R is the radius, v is the velocity, and uj is the 
angular velocity of the disk, /i is the friction coefficient, 
and the integration extends over the area of the disk with 
r vectors starting at the center. F n is the normal com- 
ponent of the force pressing the objects together at the 
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contact, in our case F n = mg, where m is the mass of the 
disk and g is the gravitational acceleration. We found it 
useful to introduce the dimensionless quantity e = v/Ruj 
with v = |v| and u = \uj\, because the friction force 
depends on v and u> only through this combination: 
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where e„ = v/u, e w = uj/oj, f = r/i?, and A\ is the 
area of the unit disk. Figure |l|(b)-(d) show local relative 
velocities on the surface of the disk for various values 
of e. Note that the local friction force does not depend 
on the absolute value of the relative velocity, only on its 
direction. After evaluating the integral in Eq. one 
gets F = —fiF n J r (e)e v , where 
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Here K(e) and E(e) are the complete elliptic integral 
functions of the first and the second kind, respectively 
p3[ . This calculation and the others below were per- 
formed using the mathematics software Maple |Q . The 
two parts of T[e) are smoothly connected at e — 1, since 
lim e _>i T{e) = 8/3ir and lim £ _^i T'{e) — 4/3tt from both 
the left and the right hand side. Here prime denotes dif- 
ferentiation w.r.t. e. The limiting values are F(0) = 
and lim e _ (00 F{e) = 1. 
The friction torque is 
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and after calculating the 
—fxF a RT(s)e u , where 
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The two parts of this function are also smoothly con- 
nected, as lim^i Tie) = 8/97T and lim^i T'(s) = 
— 4/37T from both the left and the right hand side. The 
limiting values are T(0) = 2/3 and lim e ^ 00 T(e) = 0. 
Figure || shows !F{e) and T(e), and also the T(!F) func- 
tion. This latter exists and is invertable because both 
J-(s) and Tie) are strictly monotonic functions. 

Now let us calculate how a sliding and spinning disk is 
slowing down. Assuming that only gravity and friction 
forces are acting, the scalar equations of motion are 

dv , . . . 

m— = —fimgF(e), (2) 
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FIG. 2: (a) The dimensionless friction force and torque, T 
and T, as functions of the dimensionless velocity parameter 
e. (b) The friction force and torque are coupled: the curve 
shows the possible (T, T) pairs. 



By introducing dimensionless velocities and time as v* 



v/y/Rgii, lu* — lo^J R/ g / n, and t* - 
and (j|) reduce to 

dv* 

IF = -^ (£) ' 

dt* y ' 



ty/g/R, Eqs. 
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with e — v* /lu* . As F(e) and T(e) are positive for e > 0, 
the translational and angular velocities are strictly mono- 
tonically decreasing in time, as expected. Now the ques- 
tion arises: Is it possible that any of them reaches zero 
before the other, i.e., may it happen that an initially slid- 
ing and spinning disk after some time is only sliding or 
spinning? We solved Eqs. (0) and (^|) numerically with 
many different initial conditions. The results indicate, as 
can be seen in Fig. ||(a), that v* and uj* always reach zero 
together, meaning that the disk always stops its sliding 
and spinning motion at the same moment. Before we 
show its proof, let us try to explain qualitatively what is 
happening. If the velocity is much higher than the an- 
gular velocity (v ^> Rui, i.e., e ^> 1), then the friction 
torque is negligible compared to the force, see Fig. ||(a). 
Therefore, the velocity decreases with a higher rate than 
the angular velocity, and e decreases. On the other hand, 
if the angular velocity is much higher than the velocity 
(e <C 1), then the friction torque is higher than the force, 
and e increases. Thus a negative feedback effectively 
equilibrates the sliding and spinning motion. Indeed, the 
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FIG. 3: Numerical results: (a) u*-v* trajectories with dif- 
ferent initial conditions, (b) e-v* trajectories with the same 
initial conditions as in (a). 
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FIG. 4: Function f(e) = e - J r (e)/2T(e), the r.h.s. of differ- 
ential equation ([]). It has zero value at e — and eo ~ 0.653, 
and is positive for < e < Eo, negative for eo < e. 

numerical results show this behavior, as e always tends 
to the same value, Sq « 0.653, when the motion stops 
[Fig. |(b)] . This means that lo* and v* not only reach 
zero simultaneously, but also that their ratio approaches 
a universal value, irrespective of the initial conditions. 

In order to prove that e always has this value at the 
end of the motion, we derive an autonomous differential 
equation for e from equations (Q) and (|^) using the vari- 
able transformation x = —hxu>*: 

Note that ui* — > 0, the condition of stopping, now cor- 
responds to x — > oo (with the exception of pure sliding 
motion) . For small e the right hand side of Eq. (^) van- 
ishes like /(e) rs e/4, while it behaves asymptotically for 
£ — ► oo like /(e) « — e. In between, at eo, it changes sign 
(Fig. ||). Therefore Eq. (0) has three fixed points: Two 
of them, e = and e = oo, are trivial and correspond to 
pure spinning or pure sliding motion, respectively. For 
all other initial conditions (0 < e < oo), corresponding 
to initial sliding and spinning, eo is the attractive fixed 
point, meaning that e has this value just before the disk 
stops its motion, which is what we wanted to prove. 

We also performed a simple experiment to measure the 
friction force and torque acting on a sliding and spinning 
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FIG. 5: Experimental data: the instantaneous dimensionless 
friction force (crosses) and torque (circles), acting on a sliding 
and spinning CD disk, as functions of e. One pair of data 
points is presented for each throw, taken shortly after setting 
the disk into motion. The corresponding theoretical curves 
are also displayed: T{e) (solid line) and T(e) (dashed line). 

disk. We set a standard writable CD disk (R = 6 cm), 
with its data carrying side down, into motion manually 
on a horizontal polyamid fabric surface several times, and 
recorded its motion with a Sony DCR-VX2000E PAL dig- 
ital video camera (25 images/second). Then we processed 
the images to obtain the position and the orientation 
of the disk as a function of time. We fitted the posi- 
tion and angular data with second degree polynomials to 
get the instantaneous translational and angular velocity 
and acceleration. Then, assuming that only gravitational 
and frictional forces were acting on the disk, using Eqs. 
(0) and (||), and having only the friction coefficient as a 
fit parameter, we were able to plot functions J-(e) and 
T(e), see Fig. [|. We had, however, a minor complica- 
tion: We observed that the friction coefficient slightly 
increased linearly with the number of throws, n. We 
found that /i = 0.202 + 0.00053n was a good fit, and the 
data presented in Fig. |^ was obtained using this "time- 
dependent" friction coefficient in the data processing. In- 
vestigating the disk we concluded that the reason for the 
increase of \i was probably that the film on the CD was 
gradually removed. The disk also had a 1.5 cm diameter 
hole at its center, but a numerical calculation of T(e) and 
T(e) for this geometry showed that deviation from the 
full disk case is not significant, it is much smaller than 
the accuracy of the measured data. Our experimental 
data also showed that the disk always stopped its sliding 
and spinning motion at the same moment (within error). 

We used the sliding and spinning disk as a simple, il- 
lustrative example to show how friction force and torque 
are coupled. In this case we were able to derive all re- 
sults analytically, because the local pressure is every- 
where the same in the contact area. However, in general 
the pressure distribution over the contact area will be 



4 



non-uniform. As an example, if we replace the flat disk 
by a cylinder standing on one of its flat faces, then the 
friction force leads to a torque with respect to the cen- 
ter of mass. Provided that the cylinder does not topple, 
this torque must be compensated by a pressure increase 
at the front and a pressure decrease at the rear part of 
the contact area. Therefore the spinning motion induces 
a friction component perpendicular to the translational 
motion, in the direction of v x ui. Hence, in contrast to 
the straight sliding of a flat spinning disk, the path of 
the cylinder will be curved in this direction. This resem- 
bles the Magnus effect, although the physical origin is 
completely different. 

The pressure distribution can also depend on the 
shape and elastic properties of the sliding body. For in- 
stance, if it is a sphere, linear elasticity theory predicts 
a y/l — r 2 /R 2 shaped radial pressure function Wc 
calculated the T{e) and T(e) curves numerically for this 
case and found that their qualitative behavior remains 
the same [jifl. Therefore coupling between the friction 
force and torque is still present: For large e torsion fric- 
tion is suppressed by sliding, for small e sliding friction 
gets reduced by spinning. This may explain, why the 
translational motion of a fast spinning top is hardly de- 
celerated. 

Now let us consider what impact the coupling between 
sliding and spinning motion may have on the physics of 
granular media. In one limiting case, when the particles 
are spherical and very hard, torsion friction typically can 
be neglected, as the contact area is very small, hence 
£ ^> 1. However, on the one hand, real particles are usu- 
ally non-spherical, or they may be very soft, so that the 
size of the contact area can be comparable to that of the 
particle. On the other hand, even in the case of hard 
spherical particles with contact radius r much smaller 
than the particle radius R, the coupling between sliding 
and spinning motion can have subtle consequences. As 
the torsion friction is very small in this case, one can ex- 
pect that the spinning degree of freedom can be easily 
excited. Typical sliding velocities will be comparable to 
oj n R, so that e as R/r 3> 1. In this case, the sliding fric- 
tion is basically /iF n . However, when a granular packing 
relaxes into a static configuration, the coupling between 
sliding and spinning becomes important. For example, if 
the initial sliding velocity is zero, but the spinning de- 
gree of freedom is excited (i.e., e = 0), an arbitrarily 
small force can induce sliding. 

An important extension of this work will be to inves- 
tigate the coupling between the static friction force and 
torque. As static friction is different from sliding friction, 
and its theory is somewhat more difficult, we cannot ex- 
pect that the maximum static friction forces and torques 
will lie on the friction-torque curve in Fig. |^(b). How- 
ever, we have preliminary evidence that such coupling is 



present also in the static case. We expect that the thresh- 
old torques and forces needed to turn a sticking contact 
into a sliding and/or spinning one form a curve which lies 
above the one in Fig. ||(b). In particular, this implies that 
the application of a torque at a sticking contact makes it 
easier to excite the sliding degree of freedom. For hard 
spherical spheres a very small torque will already have 
this effect. 

Finally, the fact that T{e) and T(e) depend on the 
pressure distribution in the contact area raises the ques- 
tion, whether the "inverse problem" has a unique solu- 
tion, i.e., are the experimentally accessible functions !F(e) 
and T{e) a fingerprint of the pressure distribution in the 
contact area? This would also be an interesting question 
to study in the future. 
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